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Abstract. We apply the ADM approach to obtain a Hamiltonian de- 
scription of the Einstein-Hilbert action. In doing so we add four new 
ingredients: (i) We eliminate the diffeomorphism constraints, (ii) We 
replace the densities yfg by a function ip(x,gij) with the help of a fixed 
metric \ such that the Lagrangian and hence the Hamiltonian are func- 
tions, (iii) We consider the Lagrangian to be defined in a fiber bundle 
with base space <So and fibers F(x) which can be treated as Lorentzian 
manifolds equipped with the Wheeler-DeWitt metric. It turns out that 
the fibers arc globally hyperbolic, (iv) The Hamiltonian operator H is 
a normally hyperbolic operator in the bundle acting only in the fibers 
and the Wheeler-DeWitt equation Hu = is a hyperbolic equation in 
the bundle. Since the corresponding Cauchy problem can be solved 
for arbitrary smooth data with compact support, wc then apply the 
standard techniques of QFT which can be naturally modified to work 
in the bundle. 
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1. Introduction 

QFT has been very successful in quantizing non-gravitational fields while 
ignoring the interaction with gravity. In the attempts to quantize gravity, 
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on the other hand, canonical quantization has mostly been utilized which re- 
quires to switch from the Lagrangian to the Hamiltonian viewpoint. However, 
the Lagrangian is degenerate resulting in two constraints, the diffcomorphism 
constraint and the Hamiltonian constraint. The diffcomorphism constraint 
is usually ignored since nobody knows how to handle it. The Hamiltonian 
constraint leads to the Wheeler-DeWitt equation which could only be solved 
by assuming a high degree of symmetry. 

In this paper we use canonical quantization to obtain a setting in which 
the standard techniques of QFT can be applied to achieve quantization of the 
gravitational field, i.e., gravity can be treated like a non-gravitational field. 

In order to make this approach work four new ideas had to be introduced 
in the process of canonical quantization: 

(i) We eliminated the diffcomorphism constraint by proving that it suffices 
to consider metrics that split according to 

(1.1) ds 2 = -w 2 (dx ) 2 + g ij dx i dx j 

after introducing a global time function a; . The underlying spacetime N = 
N n+1 can be considered to be a topological product 

(1.2) N = IxS 

where / C R is an open interval, So a Cauchy hypersurface, fixed for all 
metrics under consideration, and gij = gij(x ,x), x 6 So, a Riemannian 
metric. 

(ii) The volume element yjg, g = det(<?jj), is a density and it appears ex- 
plicitly in the Lagrangian and in the Hamiltonian. However, the Hamiltonian 
has to be an invariant, i.e., a function and not a density. To overcome this 
difficulty we fixed a metric x S T ' 2 (Sq) and defined the function ip by 

(1.3) = 

such that ip = (fix, gij) and 

(1-4) Vs = vVx- 

The density «/% will be later ignored when performing the Legendre trans- 
formation in accordance with Mackey's advice to only use rectangular coor- 
dinates in canonical quantization, cf. [8, p. 94]. 

(iii) After the Legendre transformation the momenta depend on x € So ■ To 
overcome this difficulty we consider a fiber bundle with base space So where 
the fibers are the positive definite metrics gij{x) over x, i.e., a fiber F(x) is an 
open, convex cone in a finite dimensional vector space. We treat this cone as 
a manifold endowing it with the DeWitt metric which is Lorentzian. It turns 
out that F(x) is globally hyperbolic. Let us call the bundle E. Each fiber 
has a Cauchy hypersurface M(x) and we denote the corresponding bundle 
by E. 
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DeWitt [3] had already the idea to consider F(x) as a Lorentzian space 
but he did not consider it to be a fiber and his DeWitt metric was no real 
tensor since it contained a density. 

The introduction of the bundle E simplifies the mathematical model af- 
ter canonical quantization dramatically. The Hamiltonian operator H is a 
normally hyperbolic differential operator acting only in the fibers which are 
globally hyperbolic spacetimes and the Wheeler-DeWitt equation is the hy- 
perbolic equation 

(1.5) Hu = 0, 

where u is defined in E. 
The Cauchy problem 

Hu = f 

(1.6) u hl =u 

D v u\ M = m 

is uniquely solvable in E with u £ C°°(E, K), K = R V IK = C, for arbitrary 

(1.7) u , Ul £C™(E,K) A f£C™{E,K). 

(iv) In view of (1.6) the standard techniques of QFT, slightly modified 
to accept the present setting, can be applied to construct a quantum field 
$g which maps functions u £ C^°(E,~Bl) to self-adjoint operators in the 
symmetric Fock space created from the Hilbert space 

(1.8) H & = L 2 (EX)- 

The quantum field also satisfies the Wheeler-DeWitt equation in the distri- 
butional sense. 

This was a summary of the contents of Sections 3-6. In the last section we 
consider the interaction of gravity with a scalar field y 

(1.9) y:N^M u 

where Mi is a complete Riemannian space. The corresponding Lagrangian 
is 

(1.10) L^-lg^viyfaAB-viy). 

It turns out that the combined fields can be treated just like gravity alone. 
The fibers have to be replaced by 

(1.11) F x Mi 

which are again globally hyperbolic, the Hamiltonian is a normally hyperbolic 
differential operator and the further reasoning is identical to the former one 
without a scalar field. 



1.1. Remark. The fibers F, and also F x Mi, are spacetimes with a past 
crushing singularity, a big bang, cf. Theorem 4.4 on page 16. 
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2. Definitions and notations 

The main objective of this section is to state the equations of Gaufi, 
Codazzi, and Weingarten for spacelike hypersurfaces M in a (n+ ^-dimen- 
sional Lorentzian manifold N. Geometric quantities in N will be denoted 
by (g a p), (Rap-ys), etc., and those in M by (<?ij), (Rijki), etc.. Greek indices 
range from to n and Latin from 1 to n; the summation convention is always 
used. Generic coordinate systems in N resp. M will be denoted by (x a ) 
resp. Covariant differentiation will simply be indicated by indices, only 
in case of possible ambiguity they will be preceded by a semicolon, i.e., for a 
function u in N, (u a ) will be the gradient and {u a p) the Hessian, but e.g., the 
covariant derivative of the curvature tensor will be abbreviated by R a /3jS;e- 
We also point out that 

(2.1) Ra/3iS;i = RaPi&-eXl 

with obvious generalizations to other quantities. 

Let M be a spacelike hypersurface, i.e., the induced metric is Riemannian, 
with a differentiable normal v which is timelike. 

In local coordinates, (x a ) and the geometric quantities of the spacelikc 
hypersurface M are connected through the following equations 

(2.2) x% = h ljV a 

the so-called Gaufi formula. Here, and also in the sequel, a covariant deriva- 
tive is always a full tensor, i.e. 



(2.3) xfj — x"j + r^x^x'J. 

The comma indicates ordinary partial derivatives. 

In this implicit definition the second fundamental form (hij) is taken with 
respect to v. 

The second equation is the Weingarten equation 

(2.4) v?=h$2%, 

where we remember that vf is a full tensor. 
Finally, we have the Codazzi equation 

(2.5) hij. k - h ik; j = R a p 1 sv a x l3 t x 1 x{ 
and the Gaufi equation 

(2.6) R^ki = -{h lk hji - huhjk} + R Q/3lS xfx^ xjxf . 

Now, let us assume that N is a globally hyperbolic Lorentzian manifold 
with a Cauchy surface. TV is then a topological product I x So, where / is 
an open interval, So is a Riemannian manifold, and there exists a Gaussian 
coordinate system (x a ), such that the metric in N has the form 

(2.7) ds 2 N = e^'i-dx® 2 + a l3 (x°, x)dx l dx j }, 
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where er.y is a Riemannian metric, tp a function on N, and x an abbreviation 
for the spacelike components (x l ). We also assume that the coordinate system 
is future oriented, i.e., the time coordinate x° increases on future directed 
curves. Hence, the contravariant timelike vector = (1, 0, . . . , 0) is future 
directed as is its covariant version (£ a ) = e 2 ^(— 1, 0, . . . , 0). 
Let M = graph U| s be a spacelike hypersurface 

(2.8) M = {{x°,x): x° =u(x), x G«S }, 
then the induced metric has the form 

(2.9) ftj = ^{-muj + try} 

where o~ij is evaluated at (u,x), and its inverse (g 1 -') = (ffy) -1 can be ex- 
pressed as 

(2.10) g ij = e~^{a ij + - — }, 

v v 

where (cr 4 - 7 ) = (ffy) -1 and 

/ - rtf, 

(2.11) 



u 



v 2 = 1 a 'u.a, = 1 - \Du\ 2 . 

Hence, graph u is spacelike if and only if \Du\ < 1. 

The covariant form of a normal vector of a graph looks like 

(2.12) {v a )=±v- 1 e^{l,-u i ). 
and the contravariant version is 

(2.13) (i/ a )= T «- 1 e-*(l,u < ). 
Thus, we have 

2.1. Remark. Let M be spacelike graph in a future oriented coordinate 
system. Then the contravariant future directed normal vector has the form 

(2.14) (i/*) =v- 1 e-^(l,u i ) 
and the past directed 

(2.15) (v a ) = ~v- 1 e^(l,u r ). 

In the Gaufi formula (2.2) we are free to choose the future or past directed 
normal, but we stipulate that we always use the past directed normal. Look 
at the component a — in (2.2) and obtain in view of (2.15) 

(2.16) er^v^h^ = —mj ~ r° UiUj - r^m - r^uj - r?-. 

Here, the covariant derivatives are taken with respect to the induced metric 
of M, and 

(2.17) -I% = e-*h ii , 

where (hij ) is the second fundamental form of the hypersurfaces {a; = const} . 
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An easy calculation shows 

(2.18) h iJ e-* = -^<7y-'0<ry ) 

where the dot indicates differentiation with respect to x . 

3. The Wheeler-DeWitt equation 

Let N = N n+1 be a globally hyperbolic spacetime. We consider the func- 
tional 

(3.1) J= f (R-2A), 

JN 

where R is the scalar curvature and A a cosmological constant. The integra- 
tion over N is to be understood only symbolically since we are only interested 
in the first variation of the functional, i.e., when a metric g = (gap) hi N is 
given, we are only interested in the first variation of J with respect to compact 
variations of <?, hence it suffices to integrate only over open and precompact 
subsets SI C N such that 

(3.2) J= [ (R-2A). 

Jn 

It is well known that, when the first variation of J with respect to arbitrary 
compact variations of g vanishes, the metric g satisfies the Einstein equations 
with cosmological constant A, namely, 

(3.3) G Q/3 + Ag aP = 0, 

where G a p is the Einstein tensor. 

When N endowed with a metric g is globally hyperbolic, there exists a 
global time function / e C°°(N) such that 

(3.4) ||^/|| 2 = .9^/a/^<0, 
N can be written as a topological product 

(3.5) N = IxS , J = (o,6)cE, 
where 

(3.6) 5 = /" 1 (c) ! a<c<b, 

is a Cauchy hypersurface and there exists a Gaussian coordinate system (x a ), 
< a < n, such that x = / and the metric g splits according to 

(3.7) ds 2 = -w 2 {dx ) 2 + gydx'da?, 
where (x l ), 1 < i < n, arc local coordinates of Sq and 

(3.8) g~ij = cjijix ^), x £ S , 

are Riemannian and w > is function. 

Without loss of generality we may always assume that € I and that c = 0. 
When there exists a time function and an associated Gaussian coordinate 
system such that (3.7) is valid we also say that x° splits the metric. 
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3.1. Lemma. Let (N,g) be a globally hyperbolic spacetime and let f be a 
time function that splits g. Let u> = (ui a p) be an arbitrary smooth symmetric 
tensor field with compact support and define 

(3.9) g(e)=g + eu 
for small values of e 

(3.10) |e| < e . 

// eo is small enough, the tensor fields g(e) will also be Lorentzian metrics 
that will be split by f. 

Proof. We shall only prove that the g(e) will be split by /, since the other 
claim is obvious. 

Define the conformal covariant metrics 

(3.H) g(e) = \(Df,Df)\ 2 g(e), 

then Df is a unit gradient field for each 17(e). Let So = / _1 (0) and consider 
the flow x = x(t,£) satisfying 

x = -Df, 

(312) *(0,{) = f, 

where £ € So- For fixed £ the flow is defined on a maximal time interval 
J = (ao,bo). If we can prove that J = I = (a,b) = f(N), then we would 
have proved that each metric g(e) satisfies 

(3.13) ds 2 = -(dx ) 2 + giJ dx d x j , 

where the gij are Ricmannian and depend smoothly on e, cf. the arguments 
in [4, p. 27]. 

It suffices to prove bo = b. Assume that 

(3.14) b < b, 

and let K be the support of uj. Then there exists to < bo such that 

(3.15) x{t,£)$K Vt>t Q , 
for otherwise there would exist a sequence (tk) 

(3.16) t k ^bo A x(t k ,£)€K 

contradicting the maximality of J, since there has to be a singularity" for 
the flow in bo- 
Thus, choose 

(3.17) t < ii < 6 , 
then 

(3.18) x(t 1 ,OeM(t 1 ) = r 1 (h), 

because 

(3.19) f(x(t,0)=t 



8 



CLAUS GERHARDT 



as one easily checks. 

Let y = y(t,£) be the flow corresponding to e = 0, then y covers TV by 
assumption and hence there exists £ € Sq such that 

(3.20) y(h,0=x(t 1 ,0- 
Then the integral curve 

V = -Df, 

V ' y(h,0=x(h,0, 

where the contravariant vector is now defined with the help g = g(0) 

(3.22) D a f = g a0 fp, 

would be a smooth continuation of x(t, £) past &o, a contradiction. □ 

The preceding lemma will enable us to eliminate the so-called diffeomor- 
phism constraint when switching from a Lagrangian to a Hamiltonian view 
of gravity. 

3.2. Theorem. Let (N,g) be a globally hyperbolic spacetime, f a time 
function that splits g with Cauchy hypersurface Sq. Let fl (s TV be open and 
precompact and assume that the first variation of the functional 

(3.23) J= / (R-2A) 

vanishes in g for those compact variations of g which can be expressed in the 
form 

(3.24) ds 2 = -w 2 (dx ) 2 + g %J dx l dx\ 

where (gij(x°,x)) is Riemannian, then the first variation of J in g also van- 
ishes for arbitrary compact variations. 

Proof. Let lu = ((d a f}) be an arbitrary smooth symmetric tensor with compact 
support in Q. The metrics 

(3.25) g(e)=g + euj, |e| < e 0j 

then satisfy (3.24) for small eo, in view of the preceding lemma, hence 

(3.26) SJ(g;g(0))=0. 
But the first variation is a scalar, hence 

(3.27) 8J{g-g(Q))=5J{g-u), 

where at the right-hand side we used an arbitrary coordinate system to ex- 
press the tensors. □ 
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We are now ready to look at the Hamiltonian form of the Einstein-Hilbert 
action following [1]. 

Let fi <e N be an arbitrary open, precompact set. Then we consider the 
functional 

(3.28) J = a N 1 { (R-2A), 

Jn 

where a at is a positive constant and we assume that there exists a time 
function / = x° in N with Cauchy hypersurface Sq = / _1 (0) and where, in 
view of Theorem 3.2, we only consider metrics of the form 

(3.29) ds 2 = -w 2 (dx ) 2 + g^dx'dx 3 , 

where w is an arbitrary smooth positive function and = gij(x°, x), x G Sq, 
Ricmannian metrics. Let us fix a metric g = {gap) as in (3.29), then we 
deduce from the Gaufi equation 

(3.30) R = H 2 - \A\ 2 +R-2R af} v a v es , 
where R is the scalar curvature of the slices 

(3.31) M{t) = {x° = t}, 
H the mean curvature of M(t) 

n 

(3.32) H = g i ^h ij =J2^ 

i=l 

where the principal curvatures, \A\ is denned by 

n 

(3.33) \A\ 2 = h ij h» =J2^l 



and where the second fundamental form hij of M(t) can be expressed as 
(3.34) = —\gijW 

where 



'ij 



(3.35) gij 



dt ' 

when we identify t with x°. 

The last term on the right-hand side of (3.30) can be written as 

(3.36) - 2R af3 v°'v f) = -2{H 2 - \A\ 2 ) + D a a a , 

cf. [7, equ. (4.60)]. Since the divergence term can be neglected the functional 
J is equal to 

(3.37) J = <*n 1 [ / i\ A \ 2 -H 2 + R - 2A}w^, 
where we may assume that 

(3.38) n = (a,b)xQ, f2<s=S a , 
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and 

(3.39) (a,b) m x°(N) = I. 
Using the relations (3.32)-(3.34) we finally conclude 

(3.40) J = a ; f [ {\G^' kl g %] g kl w- 2 + R- 2A}w^g, 

J a JO 

where 

(3.41) G^ kl = U<r f! -' + fi'fi' 1 } - gVg kl 



-l 



and 

(3-42) = ( 9ij ) 

The metric in (3.41) is known as the DeWitt metric and it is a Lorentzian 
metric as one can easily check. 

We are almost ready to define the corresponding Hamiltonian, we only 
need two adjustments. First, we like to replace the density 

(3.43) g = det( 9ij ) 

by 



where \ ~ (Xij ) is an arbitrary but fixed Ricmannian metric in iSo and where 
ip is now a function 

(3.45) = = ^|- 

Second, the Riemannian metrics gij(t, ■) are elements of the bundle T°' 2 (Sq). 
Denote by E the fiber bundle with base So where the fibers consists of the 
Riemannian metrics (fly). We shall consider each fiber to be a Lorentzian 
manifold equipped with the DeWitt metric. Each fiber F has dimension 

n(n + I) 

(3.46) dim F = — '- = m + 1. 

Let (£ a ), < a < to, be coordinates for a local trivialization such that 

(3.47) n 

is a local embedding. The DeWitt metric is then expressed as 

(3.48) Gab = G ij ' kl g ijia g k i, b , 

where a comma indicates partial differentiation. In the new coordinate sys- 
tem the curves 

(3.49) t^ gij (t,x) 
can be written in the form 

(3.50) 
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and we infer 

(3.51) G' J - l, 'l,,ilu = G ab ti b - 
Hence, we can express (3.40) as 

(3.52) J = [ [ a- l {\G ab ti b w' l V + (R - 2A)wip} 



where we now refrain from writing down the density explicitly, since it 
does not depend on (gij) and therefore should not be part of the Lcgendre 
transformation. Here we follow Mackey's advice in [8, p. 94] to always con- 
sider rectangular coordinates when applying canonical quantization, which 
can be rephrased that the Hamiltonian has to be a coordinate invariant, 
hence no densities are allowed. 

Denoting the Lagrangian function in (3.52) by L, we define 

or i 

(3.53) n a = -j- = <pG ab — i'w- 1 

d£ a 2aAr 

and we obtain for the Hamiltonian function H 

(3.54) = (2^Z bw ~ 1 ) waN ~ a x( R ~ 2A ^ W 
= <p~ 1 G ab iT a 7r b wa N - a N 1 (R - 2A)ipw 

= Hw, 

where G ab is the inverse metric. Since w is an arbitrary function we obtain 
the Hamiltonian constraint 

(3.55) H = a N ip- 1 G ab ir a w b - a N \R - 2A)tp = 0. 
Applying canonical quantization, by setting h = 1, we replace 

1 d 

(3.56) 7T a = TT a (x) ->• -- 



i d£ a (x)' 

where (£ a (x)) are points in the fiber over x £ So and 

< 3 - 57 ' WIS) 

denotes partial differentiation in the fiber over x. 

Each fiber can be viewed as a Lorentzian manifold equipped with the 
metric 

(3.58) a^ipGab- 

After quantization the Hamiltonian function is transformed into the hyper- 
bolic differential operator 

H = -A-a- N \R-2A)y 
(3.59 , 

= D- a N 1 {R-2A)ip 
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where □ is the d'Alembertian operator for that metric. 

We want to emphasize that H is defined in the bundle E acting only in 
the fibers. Hence, let u be a smooth function defined in E 

(3.60) u = u(x,£(x)) = u(x,gij(x)), 
then the Hamiltonian condition (3.55) takes the form 

(3.61) Uu- a N 1 (R-2A)ipu = 0. 

Since each fiber is equipped with the Lorentz metric (aJ^tpGab) there is a 
natural volume clement 



(3.62) A /|det(a^VGab)|^ 
and for functions C^°(E,K), where 

(3.63) K = C V K = R, 
we can define the scalar product 

(3.64) (u,v)e = uv, 

J s JF(x) 

where the volume element of the fiber F(x) is given in (3.62) and that of 6>o 
is defined by the metric Xij 

(3.65) y/xdx. 

ft is immediately clear that the hyperbolic differential operator is formally 
self-adjoint, and since each fiber is globally hyperbolic, which will be proved 
in Section 4, the Cauchy problem for the Wheeler-DeWitt equation 

(3.66) Hu = 

can be uniquely solved in E for given initial values on a Cauchy hypcrsurface, 
and the standard techniques of QFT can be applied to quantize the fields u 
in (3.66) which represent gravitation. 

This procedure will be outlined in the sections below. Let us also mention 
that the constructed quantum field <S> ^,E will be a Cauchy hypersurface in E, 
will satisfy the Wheeler-DeWitt equation in a distributional sense indicating 
that the QFT approach is an accordance with the Hamiltonian constraint. 



4. The fibers are globally hyperbolic 

DeWitt already analyzed the fibers in [3] , though he did not look at them 
as fibers. Some of the ideas that we shall use in the proofs of the lemma and 
the theorem below can already be found in DeWitt's paper as we discovered 
after studying the literature. 



4.1. Lemma. Let F be a fiber, then F is connected and 
(4.1) t = logy 
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is a time function satisfying 

(4.2) a N <p- 1 G ab T a T b = - ., - ^ a N ip-\ 

4(n — 1) 

Proof. F is obviously connected, since F is a convex cone in the vector space 
defined by the symmetric covariant tensors of order two. 

To prove (4.2) we use the original coordinate representation g^ and con- 
clude 

and hence 

(4.4) <Wt* - -j^. 

where 

(4-5) G ijM = ^{gik9fl + 9a9jh} - ^j9ij9u 

is the inverse of , hence the result . □ 

4.2. Theorem. Each fiber F is globally hyperbolic, the hypersurface 

(4.6) M = {tp = 1} = {t = 0} 

is a Cauchy hypersurface and in the corresponding Gaussian coordinate sys- 
tem (£ a ) the metric aJ^ipGab can be expressed as 

(4.7) ds 2 = ^ n ~ l) a~ N \{-dT 2 + G AB di A d^ B }, 
where 

(4.8) t = £° A -co < t < oo 

and (^ A ), 1 < A < m are local coordinates for M . The metric Gab *s also 
static, i.e., it does not depend on r. 

Proof, (i) Let r be as in Lemma 4.1, then t(F) = M and in the conformal 
metric 

(4.9) G 



ab = a N (fi 1 — —(a^tpGab) 

4(n — 1) 



t q is a unit gradient field in view of (4.2). 

(ii) The hypersurface M in (4.6) is therefore spacelike and has at most 
countably many connected components. 

Consider the flow 

i= -Dt = -(G ab n) 
(4.10) 5 1 ' 

£(o,C) = C, CgM. 
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It will be convenient to express the flow in the original coordinate system, 
i.e., 

4(rc-l) 

(4.11) 913 ~ n l3M ' 

9ij(0,0 = C = 9ij, 

where Gij u is the metric in (4.5). The flow exists on a maximal time interval 
J C . 

From (4.5) we obtain 

GkX 1 f i kl 
ij,kl T = ^yiiMQ 

(4-12) 1 J_ 

_2 - 9jj1 n-l>- 2{n-l) 9l] ' 

hence 

(4-13) //,, = l9ij- 

Let (77*) € T^°(Sq) be an arbitrary unit vector with respect to the metric 
Xij, then 

(4.14) (gijvWY = hh'in 1 

leading to 

(4.15) g^rfrf = g^rfrfe^, 

thus the eigenvalues of with respect to Xij are uniformly bounded from 
above and strictly bounded against zero when |£| < const. Moreover, 

(4.16) T( 9ij ) = t 
from which we conclude 

(4.17) J ( = R. 

If M would be connected, then we would have proved that F is product 

(4.18) F = M x M 

and that the metric would split as (4.7). However, if M had more then one 
connected component, then the corresponding cylinders defined by the flow 
would be disjoint and hence F would not be connected. 

(iii) Let (£ a ), < a < to, be the corresponding Gaussian coordinate system 
such that 

(4.19) C° = t = t 

and (£ ), 1 < A < m, are local coordinates for M. Let gij(£, a ) be a local 
embedding in the new coordinate system, where the ambient metric should 
be the conformal metric up to a multiplicative constant, i.e., we consider 

(4.20) Qij.kl = l{gijgkl + gilgjk} _ gij g M 

to be the ambient metric such that 

(4.21) Gab = G ij ' kl g ij!a g k i, b . 
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The metric splits and we claim that 

(4.22) G AB = G ij ' kl g ij>A g k i,B 
is stationary 

(4.23) jGab = 0. 

To prove this equation we observe that the normal to M(t) = {t = t} is a 
multiple of g^ , cf. (4.3), hence 

(4-24) gV giitA = 

for gij(t, £ ) is a local embedding of M(t) from which we deduce 

(4.25) Gab = §{s*V' + <J ,! <J , 'l<J, J .A<nu:- 

Differentiating this equation with respect to t we infer, in view of (4.13), 

jGab = -l{g ih 9 il + g iX '■9 jk }9i j ,A9H,B 

(4-26) , 2 f„ik„jl , „il„jk\ „ 

y ' + nl9 9 +.9 9 \9ij,A9kl,B 

= 

where we also used 

(4-27) 9 ij = ~l9 ij - 

(iv) Finally, we want to prove that M — M(0) is a Cauchy hypcrsurface 
and hence F globally hyperbolic, cf. [9, Corollary 39, p. 422]. It suffices to 
prove this result for a conformal metric G a b where 

(4.28) ds 2 = -dr 2 + G AB d£, A d£ B 

and Gab is stationary. 

Gab is the metric of M. In case n = 3 DeWitt proved in [3, Remarks 
past equ. (5.15)] that M is a symmetric space and hence complete. DeWitt's 
proof in [3, Appendix A] remains valid for n > 3. We shall only use the fact 
that M is complete; in Lemma 4.3 below we shall give a second proof which 
does not rely on DeWitt's result. 

Let 7(s) = (7°(s)), s £ /, be an incxtcndiblc future directed causal curve 
in F and assume that 7 does not intersect M. We shall show that this will 
lead to a contradiction. It is also obvious that 7 can meet M at most once. 

Assume that there exists sq £ I such that 

(4.29) r( 7 (s )) < 

and assume from now on that sq is the left endpoint of I. Since r is continuous 
the whole curve 7 must be contained in the past of AI . 
7 is causal, i.e., 

(4.30) G AB i A i B < |7°| 2 
and thus 



(4.31) VG AB i A i B < 7° 
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since 7 is future directed. Let 

(4.32) 7 = (7 A ) 

be the projection of 7 to M, then the length of 7 is bounded 

(4.33) £(7) = J VG AB i A i B < J 7° < -7°(*o). 

Hence, 7 stays in a compact set since M is complete and the timclikc coeffi- 
cient is also bounded 

(4.34) 7°(s ) < 7 °(s) < Vs G J, 

which is a contradiction since 7 should be inextendible but stays in a compact 
set of F. □ 

4.3. Lemma. The hyper surf ace M — Mix) is a Cauchy hypersurface in 
F(x). 

Proof. As in the proof above we consider an inextendible causal curve 7 and 
look at the projection 7 given in (4.32) which has finite length, cf. (4.33). 
Then, it suffices to prove that 7 stays in a compact subset of M. 

Representing 7 = 7(5), s 6 I = [so,b), in the original coordinate system 

(9ij) 

(4.35) j=( gij (x,s)) = ( gij (s)) 
we use (4.25) to deduce 

(4-36) L( 7 ) = < - 7 °( So ), 

where 

(4-37) ;/,, " 'r ^ii.jilu. 

from which we infer, in view of [6, Lemma 14.2], that the metrics (gij(s)) are 
all uniformly equivalent in I and converge to a positive definite metric when 
s — > b. Hence, the limit metric belongs to M and 7 stays in a compact subset 
of M. □ 

4.4. Theorem. Each fiber F is a spacetime with a past crushing singu- 
larity, a big bang. It is past-timelike incomplete and the volume of the slices 
M(t) converges to as r — > — 00, or equivalently, when tp — > 0. 

Proof. The mean curvature H of the slices M{t) with respect to the past 
directed normal is 



(4.38) i? = -^-^V^- 1/2 , 

as one easily checks. Moreover, let v be their past directed normal vector 
field, then 

(4.39) Rabv a v h = 0, 
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where (R a b) is the Ricci tensor of the ambient space, cf. [3, equ. (A35) 
in Appendix A]. Hence, Hawking's singularity theorem can be applied to 
conclude that F is past-timelike incomplete. Furthermore, the sectional cur- 
vatures squared of the Riemann tensor tend to infinity if the singularity is 
approached, see [3, Appendix A]. □ 

5. Solving the Wheeler-DeWitt equation 

In this section we want to solve the Wheeler-DeWitt equation (3.62) on 
page 12 in the bundle E. Let us first recall some well-known definitions 
and results for hyperbolic differential operators of second order in globally 
hyperbolic spacetimes. 

5.1. Definition. Let N be a smooth, globally hyperbolic spacetime. A 
linear differential operator P of order two in N is said to be normally hyper- 
bolic if it can be expressed in the form 

(5.1) P = -A + b a D a + c, 

where A is the Laplace operator of the underlying metric, (b a ) a smooth 
vector field and c a smooth function. 

5.2. Theorem. Let N be globally hyperbolic, M (Z N a Cauchy hypersur- 
face with future directed normal v , P a normally hyperbolic operator, uq, u\ 
resp. f functions in C£°(M,K) resp. C^°(N, K) ; then the Cauchy problem 

Pu=f, 

(5.2) u\ M =u , 

u a v a \ M = «i, 

has a unique solution u G C°°(N, K) such that 

(5.3) suppu C J N (K), 
where 

(5.4) K = suppuo U suppui U supp/ 
and 

(5.5) J N {K) = J^(K) U J N (K), 

these are the points that can be reached by causal curves starting in K . More- 
over, u depends continuously on the data {uQ,ui,f) with corresponding es- 
timates, namely, for any compact sets K, K\ C N and Kq C M and any 
m G N there exists ml G N and a constant c = c(m, m', K, Kq, K{) such that 

(5-6) \u\ m ,K < c{\u \ m >, Ko + |lii| m /,jr + l/lm'.ifi}, 

where u is a solution of the Cauchy problem and uq, u\ and f have support 
in the respective sets Kq and K\ . 



A proof is given in [2, Theorem 3.2.11, Theorem 3.2.12]. 
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5.3. Remark. A corresponding result, when the data have no compact 
support, is also valid, cf. [5, Corollary 5 on p. 78]. In this case the estimates 
are of the form that for any compact K C N there is a compact K' C N 
such that for any m E N there exists m! € N and a constant c such that 

(5.7) \u\m,K < c{\u \ m i.K>nM + \ui\ m ',K'nM + \f\m',K'} 

for all solutions u of the Cauchy problem with smooth data (uq, Ux, /)• 



The Hamilton operator in (3.59) on page 11 is certainly normally hyper- 
bolic in each fiber, hence, given 

(5.8) uo,*,! G C™(E,K) 
and 

(5.9) feC™(E,K) 
the Cauchy problem 



(5.10) 



u a v 



Hu 


= f 


u \m 


= u 






■ a \ 

\M 


= Ul 



has a unique solution 

(5.11) u = u(x, £(x)) 

where x £ So, £ a (x) are local coordinates for the fibers F(x) and E is the 
bundle with base space So and fibers M(x), such that 

(5.12) u(x,-) € C°°(F(x),K). 

We shall prove in the theorem below that the solutions are also smooth in x 
such that u g C°°(E,K). 



5.4. Theorem. The solution u of the Cauchy problem (5.10) with the data 
given in (5.8) and (5.9) belongs to the space C°°(E, K) and satisfies 

(5.13) supp«cJ £ (X)= |J J F{x) {K{x)), 

xdSo 

where 

(5.14) K = suppito U suppui U supp/ 
and 

(5.15) K{x) = Knir' 1 ^), x e S , 
and it is the projection from E to So- 
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Proof. We may work in a local trivialization (£, U) of E 

(5.16) £ : t/ x (Q c R m+1 ) -> 7r _1 (l7) 
such that 

(5.17) € = e(ar,C), C G fi. 

Since the fibers are manifolds £ can be expressed in coordinates £ = (£ a (x, £))■ 
A function u = can then also be written in the form 

(5.18) u(a;,0 0)- 

The coordinates £ = (C a ) are then coordinates for the fibers F(x), x E U, 
such that the differential operator P has the form 

(5.19) Pu = -G ab u ab + b a u a + cu, 

where the derivatives are partial derivatives with respect to £ and the coef- 
ficients and it depend on (x,C). We shall prove that u(x, Q is smooth in x 
by using the difference quotient method. For simplicity we use the symbol £ 
instead of £ to denote the coordinates. 

(i) We shall first prove that u is Lipschitz continuous in x. Let x = (x l ), 
fix 1 < k < n and let h ^ be a real number. By a slight abuse of notation 
we define the n-tuple 

(5.20) h= (0,...,0,M,...,0), 
where only the fc-th component is different from 0. 

Let 

(5.21) v = h^iuix + h, ■) - u(x,-)), 
then v solves the differential equation 

(5.22) Pv = f h , 

where fh(x, •) € C°°(F(x),K) such that for any given compact K C E and 
p 6 N 

(5.23) |//i|p,irnF (x) < c = c(p, i^) V < e 
uniformly in x, with initial values 

(5.24) h- 1 {u (x + h,-)-u(x, ■)) 
and 

(5.25) ft _1 (ui(a; + ft, •) - u x {x, ■)). 

To verify (5.22), (5.23) we only consider the main part of the differential 
operator 

(5.26) G ab (x + h)u ab (x + h)- G ab (x)u ab (x), 
where u ab are partial derivatives and use the algebraic identity 

G ab {x + h)u ab {x + h)- G ab {x)u ab {x) 

= G ab (x){u(x + h)- u(x)} ab + {G ab (x + h)- G ab (x)}u ab (x + h). 
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Similar expressions occur when we look at the lower order terms. Hence 
(5.22) and (5.23) are proved. 

Then we apply the a priori estimates (5.7) to conclude 

(5.28) \v\ p ,KnF(x) < c(p,K) 

for any compact subset of K C E and any p £ N independent of h and x G U. 
Therefore u and its derivatives with respect to £ are Lipschitz continuous in 
x. 

(ii) Next we shall prove that u is of class C 1 in x. Differentiate the equation 

(5.29) Pu = f 

with respect x k , 1 < k < n, and pretend that u can be differentiated, then 
we obtain 

(5.30) Pu k = /, 
where 

(5.31) * - £ 

and f(x, •) € C°°(.F(x), K) and Lipschitz continuous in x. 
Let w be a solution of the Cauchy problem 

(5.32) Pw = f 
with initial values 

(5 33) duo A dui_ 

dx k dx k ' 

w = w(x,£) is of class C°° in ^ and Lipschitz continuous in x because of the 
arguments in (i). When we look at the difference 

(5.34) P(v -w) = f h - f, 

where v is defined in (5.21), we deduce that the right-hand side converges to 
zero in the limit h — > 0. 

The same result is also valid for the difference of the initial values, hence 

(5-35) g^=w 

and u is of class C 1 in x as well as its derivatives with respect to £. 

(iii) We can now prove by induction that u is in C°°(E, K). □ 



5.5. Corollary. Let P be as in the preceding theorem. For (x,£) G E let 
F± be the fundamental solutions with respect to P with past-compact resp. 
future- compact support at (x,^) and let G± be the advanced resp. retarded 
Green's operators of P. Then, for any u € C^°(i?,K), the mappings 

(5.36) (x,0 ^ F±(x,0[u} 
are smooth and 

(5.37) G±(u) E C°°(E,K). 
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Proof. ,,(5.36)" We only prove the claim for F + . There exists a a function 
x(it), which is given as the solution of the Cauchy problem 

(5.38) P*x{u) = u 

with vanishing initial values on a Cauchy hypersurface, where P* is the ad- 
joint of P — we do not assume P to be formally self-adjoint — such that 

(5.39) F + (s,OH=x(«), 

cf. [5, equ. (3.10) on p. 79]. Since P* is also normally hyperbolic and smooth, 
X( u ) £ C°°(E, K)) because of the preceding theorem and the claim is proved. 

,,(5.37)" The Green's operators are defined by 

(5.40) (GzpuXz.O = F±(x,0[v] V(z,0 £ E, 

hence the result because of (5.36). □ 

6. Applying QFT techniques to the Wheeler-DeWitt equation 

The main difference between the standard QFT setting and the present 
situation is that we consider a bundle where the fibers are globally hyperbolic 
spacetimes and the hyperbolic operator is acting in the fibers. However, 
because of the results in Section 5, we can solve the Cauchy problem for the 
Hamiltonian operator H in the bundle E, hence there exist the advanced and 
retarded Green distributions G+ and G_ for H such that 

(6.1) G± : C™(E,K) -> C°°(E,K) 

(6-2) H o G± = G± o H | C?W) = id Cr{E , K) 

(6.3) supp(G±M) c Jf(suppit) Vwe C™{E,K) 
and 

(6.4) supp(G_u) C J B (suppu) Vw € C^°(E, K). 
H is formally self-adjoint 



(6.5) 



(Hu, v)e — I I Huv 
/So Jf 



/ / uHv = (u, Hv)i 

JSa JF 



/S JF 

where the volume element in F is given by the Lorentzian metric and that 
in So is defined by the previously used metric x- 

There are two ways to construct a quantum field given a formally self- 
adjoint normally hyperbolic operator in a globally hyperbolic spacetime or, 
in our case, in the bundle E. One possibility is to define a symplectic vector 
space 

(6.6) V = C™(E,K)/N(G), 
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where G is the Green's distribution 

(6.7) G = G+-G-. 
Since 

(6.8) G* = -G 
the bilinear form 

(6.9) / / (Gu, v) u, v S V 

JSo J F 

is skew-symmetric, non-degenerate by definition and hence symplectic, and 
then there is a canonical way to construct a C*-algcbra. 

The second method is to use a Cauchy hypersurface to define a quantum 
field in Fock space. We shall show that this is also possible in our case. 

First let us prove the following lemma: 

6.1. Lemma. For all C^°(E,K) there holds 

(6.10) / f(u,Gv)=[ [ {{D v (Gu),Gv)-{Gu,D v {Gv))}, 
J s J f Jso Jm 

where v is the future normal to M and the scalar product is the standard 
scalar product in K. 

Proof. Let F = F(x) be an arbitrary fiber, then 

(6.11) [{u,Gv)=[ (u,Gv)+ [ (u,Gv), 

J F J F + J F- 

where 

(6.12) F + = {£° > 0} A F_ = {£° < 0}. 
Now, in F + we have 

(6.13) HG-u = u 
and 

(6.14) HGv = 0, 

in view of (6.2). Moreover, because of (6.4), we have 

(6.15) supp (G-ii) n F + is compact, 
since 

(6.16) supp(G_u)nM is compact. 
Hence, we obtain by partial integration 

(6.17) f (HG-u,Gv) = - f {D v G-u,Gv)+ [ (G-U,D v Gv). 
Jf + Jm Jm 

A similar argument applies in F- by looking at 

(6.18) HG+u = u 
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leading to 

(6.19) \ (HG+u,Gv)=[ {D v G+u,Gv)- [ {G + u,D v Gv). 

JF- JM JM 

Adding these two relations implies the result in the fiber. □ 

We now define the Hilbcrt space H E which is used to construct the sym- 
metric Fock space, namely, we set 

(6.20) H E = L 2 {E,<£) 
with the standard scalar product 



(6.21) (u,v) E = / / («,«)■ 

J S JM 

We denote the symmetric Fock space by F(Hfi). Let be the correspond- 
ing Segal field. Since G* = — G we deduce from (6.16) and from a similar 
result for G+ 

(6.22) G*u\ & G C C °°(£,R) c H & . 
We can therefore define 

(6.23) <P E (u) = e(i(G*u) l& - D v {G*u) h ). 

From the proof of [2, Lemma 4.6.8] we conclude that the right-hand side 
of (6.23) is an essentially self-adjoint operator in ^(H^). We therefore call 
the map <1> E from C^°(E,R) to the set of self-adjoint operators in ^(H^) a 
quantum field for H defined by E. 

6.2. Lemma. The quantum field <P E satisfies the equation 

(6.24) H<P E = 
in the distributional sense, i.e., 

(6.25) (H<P E ,u) = (<P E ,Hu) =<P E {Hu) =0 Vti£ C^°(E, M). 

Proof. In view of (6.2) there holds 

(6.26) G*Hu = Wue C™(E,R). 

□ 

These remarks should suffice to indicate that the techniques of QFT can 
be applied in the present case. For further details we refer to [2, Chapter 
4.7]. 
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7. Gravity interacting with a scalar field 

In this section we demonstrate that other fields can be added to the func- 
tional and that the general procedure, define a fiber bundle where the fibers 
are globally hyperbolic spacetimes, can naturally be applied — at least in the 
case of a scalar field. 

Thus, let Mi be a complete Riemannian space of dimension n\ and let 

(7.1) y : N -> Mi 

be a scalar field map; in local coordinates we have 

(7.2) y = (y A (x a )). 

N is the Lorcntzian manifold we looked at in Section 3. 
We shall consider the Lagrangian 

(7-3) L^-lg^ySvggAB-Viv), 

where (gAB) is the fixed metric in Mi and V(y) a smooth potential. 

After introducing the time function x° and the Cauchy hypcrsurface iSo 
L\ can be written in the form 

(7.4) Lr = \w- 2 g ABy A y B - ^yfyfgAB - V(y). 

We still have to multiply this Lagrangian with w^fg. Introducing the function 
tp as in (3.45) on page 10, we obtain 

(7.5) Li = Iw^gAB^y 8 ^ - ^g lJ yfyf gabwip - V(y)w<p, 

where we omit yfx as before. 

Similar as in Section 3 we regard the mapping y as a section in a bundle 
with base manifold Sq and fibers Mi. Since the metrics gij are also looked 
at as the section of a fiber bundle, we define the bundle E as a fiber bundle 
with base space So and fibers 

(7.6) F(x) x Mi, 

where the fibers are a metric product. 
Defining 

(7 - 7) PA = W 

we obtain 

8 ^ PA dy~A~ ~ Ll = 55ab(2/ j4 w"V)(2/ S ^"V)w" 1 

+ \g l ' J 'vtvf ! 9ABWip + V(y)wip. 

Hence the combined Hamiltonian constraint has the form, compare with 
(3.55) on page 11, 

H = a N ip^ l G ah ^ a ^b + l<p~ 1 g AB pAPB 

+ \g lJ y?yfg A BV ~ a£{R - 2A) V + Vip = 0. 
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The common product metric looks like 

(Q 'at 1 Gab \ 

N 
2g AB J 

After canonical quantization the Hamiltonian operator has the form 

(7.11) H = -A- a N HR - 2A)<p + ^yfyfg AB <P + Vtp 
and the Wheeler-DeWitt equation is 

(7.12) Hu = 0. 
The term 

(7.13) 9 ij yfyf9AB 

is well defined since 

(7-14) (e(x),y A (x)) = ( gij ,y A (x)) 

is a smooth section in the bundle E. 
The fibers 

(7.15) F x Mi, 

equipped with the product metric (7.10), arc globally hyperbolic spacetimcs; 
the former Cauchy hypcrsurface M has to be replaced by 

(7.16) M x Mi 

and the Gaussian coordinate system associated with the new Cauchy hyper- 
surface is 

(7.17) (C,y A ). 

The proof that the fibers are globally hyperbolic is identical to the one 
given in Section 4. The previous considerations in Section 5 and Section 6 
are therefore valid in the present setting leading to a unified quantum theory. 
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